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linear) curve select the integrable dynamics of the Ablowitz-Ladik hierarchy of evolution equations: 
i) the set of points describing the discrete curve lie in the sphere 5*^; ii) the distance between any two 
subsequent points does not vary in time; iii) the dynamics does not depend explicitly on the radius 
of the sphere. These results generalize to a discrete context our previous work on continuous curves 
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1 Integrable discrete dynamics and geometry 



During the last two decades many important classes of nonlinear equations have been solved. These 
integrable nonlinear equations appear in very different forms: they are ordinary differential equations 
(ODE's) (autonomous and not)||], partial differential equations (PDE's) 1^, algebraic and functional 
equations [^ . It was also possible to construct integrable discrete analogues of all the important nonlinear 
ODE's and PDE's il [i il [§ i) [§ ; in particular, Ablowitz and Ladik (AL) have constructed an 
integrable discrete analogue of the AKNS hierarchy ||l^ , associated with the following discrete eigenvalue 
problem 



[ ^n. 1 




. ^n. . 





c 



qk 






(1) 



where the complex fields qk and r^ depend on the discrete variable fc (fc G Z), and the complex parameter 
C plays the role of spectral parameter. 

The interesting connections between the differential geometry of submanifolds and integrable nonlinear 
PDE's, which have been shown many times throught the years (the literature on this subject is very large 
and we refer to our previous paper |ll| for a historical review), can also be found at a discrete level. For 
instance, discrete pseudospherical surfaces and the discrete analogues of constant mean curvature surfaces 
are described by integrable discrete analogues of the sine-Gordon (SG) |12 and sinh-Gordon equations 
p3[ . Such discretizations were found by adapting the main geometric properties of the continuous surfaces 
to a discrete level. 

In a recent paper ||l| we have shown that the following three elementary geometric properties of the 
motion of a curve select, among all possible dynamics, integrable dynamics: 

Property 1. The motion of the curve takes place in the A'^-dimensional sphere of radius R, denoted by 

5^(i?), iV > 1. 

Property 2. The curve does not stretch during the motion. 

Property 3. The dynamics of the curve does not depend explicitly on the radius R. 

These dynamics are described by integrable commuting nonlinear evolution equations in 1+1 (one spatial 
and one temporal) dimensions for the geodesic curvatures of the curve. More precisely, for N — 2 we 
have obtained the modified Korteweg - de Vries (niKdV) hierarchy for the curvature n. For iV = 3 
we have obtained a hierarchy of evolution equations for the geodesic curvature k = ki and torsion 
T ~ K2, which can be transformed, through the Hasimoto transformation, into the nonlinear Schrodinger 
(NLS) hierarchy. For iV > 3 we have obtained integrable multicomponcnt generalizations of the above 
hierarchies. Also integrable evolution equations "with sources" (with singular dispersion relation) have 
been selected by Properties 1-3. 

The goal of this paper is to construct large classes of geometrically meaningful discrete evolution 
equations applying the above Properties 1-3 to the case of a discrete curve, where by discrete curve we 
mean just a sequence of points. The obvious discrete generalization of Properties 1-3 is the following. 

Property i'. The motion of the discrete curve takes place in the iV-dimensional sphere S^{R) (N > 1) 
of radius R. 

Property 2\ The distance between any two subsequent points does not vary in time. 

Property 3\ The dynamics of the discrete curve does not depend explicitly on the radius R. 

We will consider the cases N = 2 and iV = 3, showing that Properties 1 '-5' generate the AL hierarchy of 
evolution equations. In particular we will prove that Property f for N = 2 and iV = 3 is equivalent to 
the spectral problem (|l|) in the reductions r^ = —qk G R and r^ — —qk respectively. 

In this paper we concentrate our attention on continuous evolutions, i.e. time i is a continuous 
variable; we postpone to a subsequent work the study of integrable discrete dynamics of a curve. 



2 The motion of a discrete curve on S'^{R) 

The problem of a curve on S'^ is a special of case of that of a curve in S^, but because of its pedagogical 
importance we have decided to treat it first. 

2.1 The discrete curve on S'^{R) 

Generalizing the problem of a discrete curve on the plane (see, eg. [^), we derive now the set of equations 
which describe a discrete curve on S^{R). We indicate by r^, (fc 6 Z) the position vector of the k^^ point 
of the sequence pointing from the centre of the sphere. The unit vector in the direction of f^ is denoted 
by ffc. The oriented plane tt^ is given by the ordered pair (ffe , rfc+i). By t^ (tj^,) we denote the unit 
vector of the plane tt^ orthogonal to ffc (ffc+i) with the direction given by the orientation of the plane TTfc. 
Both tfe and t'/^ are tangent (in the points r^ and fk+i correspondingly) to the big circle passing through 
the points r^ and r^+i. 

In the following we assume that the segments of geodesies which join subsequent points have the same 
length A for each pair of the sequence. The angle a is one-half of the angle AA (where A = l/R ) between 
the vectors r^ and Vk+i'- 

a = A/{2R) . (2) 



Together with the basis (f fc, tfe) of the plane tt^ we will use also the new basis (fj^,, a^) obtained from the 
old one by rotation of the angle a 

r'^ = ffccosa - tfcsina , 

afc — tfccosa + ffcsina 



We remark that 



rfc+i - rfc 
afc = -p ^ . (4) 



We define the normal vector n^ by fife = f ^ x tfc, where the cross-product x is consistent with the given 
orientation of the ambient R'^. The three vectors (ffe, tfc, fife) form the orthonormal basis of R^ in the 
point Ffe of the discrete curve. 

The oriented tangent plane Tfc is given by the ordered pair of vectors (tfc, fife). Since ffe _L Tfc and 
ffc_Lt'^_-^, then tjj._.|^ G Tfc. By i^fc_i (|(/3fe_i| < it) we denote the angle between t'f._^ and tfc measured 
according to the orientation of the plane Tfc; ipk-i is the angle of geodesic curvature. We remark that 
ipk-i is also the angle between the planes nk-i and TTfc. 



In the continuous limit: A — > and Ak — > s, where s is the arc-length parameter of the corresponding 
curve, the angle of geodesic curvature reduces to the familiar geodesic curvature k{s) in the following 
way: 

— — > k{s) . (5) 

The transition from the basis (ffc, tfc, fifc) to the basis (ffc+i, tfc+i, fifc+i) can be obtained by the superpo- 
sition of two rotations: 

i) rotation in the plane TTfe of the angle AA = 2a 

(ffc, tfe, fifc) — > (ffc+i, tfc, fifc) 



ii) rotation in the plane Tk+i of the angle ipk 

(f/c+i,tfe,fife) — > (ffc+i,tfc+i,nfe+i) 
The full rotation can be written in matrix form as: 



tfe+i 
n/c+i 



cos(AA) sin(AA) 

- cos ipk sin(AA) cos ipk cos(AA) sin ip^. 
sinipfc sin(AA) — sin(y9fc cos(AA) cos^s^ 



rfc 
tfe 
fifc 



(6) 



2.2 The spinor representation of the hnear problem 

To obtain the 2x2 version of the linear problem ^ we use the standard representation of rotations in 
R'^ in terms of S'{/(2)-valued matrices [^. From this point of view R'^ can be considered as the space 
spanned by the matrices ei — 102-, £2 — — icri, £3 = ia^, (here ai denote the standard Pauli matrices): 



vGR3 <=» cr(v)=^v'£/ ; ||v|p = det cr(v) 



(7) 



i=i 



where v' are the Cartesian coordinates of the vector v. We recall that any rotation O in R'^ can be 
represented by a matrix S S SU{2) by the formula 

w = Ov -^ cr(w) = S*" V(v)S' . (8) 

Suppose we are given a sequence of orthonormal frames Tk = (e/j, ffc, gfe), k G Z, and the corresponding 
matrices iS^ 



■'k — ^k ^^k J 

where £ = (£i,£2,£3)- When matrices Ak are defined by the relation 

Sk+i = AkSk 
then they satisfy also 



SkJ-k+iSi^ — Af, £Ak ■ 



For example, the matrix 



Ak 



e'" 
e~^° 



represents the rotation in the plane (e^;, ffc) of the angle 2a, and the matrix 

Ak-- 



cos((^/2) sin((/3/2) 
— sin{(p/2) cos{(p/2) 



(9) 
(10) 

(11) 
(12) 

(13) 



represents the rotation in the plane (fj,, gj.) of the angle ip. 

When the sequence of orthonormal frames is the sequence of the Frenet frames in the points of the 
discrete curve: J-k = (f j., t^, n^;), then the corresponding spin matrices Sk are subjected to the relation 



SkA 



e*"cos(</7fe/2) e-'"sm{ipk/2) 
-e'"sin(</7fc/2) e~'°cos(</jfe/2) 



1 



VT+^ 



<qk C' 



Sk 

Sk —'■ AkOk 



(14) 



with 

gfe-tan(^fe/2) , C = e" . 

We remark that this spectral problem is equivalent to the AL spectral problem 



Sk+i — 



C qk 

-qk C 



Sk 



(15) 



(16) 



in the reduction rk = —qk G R- To show this we relax the (unessential) condition det Sk = 1 to remove 
the cos((/3fc/2) diagonal term and apply a gauge transformation which corresponds to the change of the 
basis of the plane tt^ from (ffe,tfe) to (f'j.,afc): Tk = (ffc,tfe,nfe) ^ Tk = (f'j.,afe,nfe). Then 



where 



and 



Sk = Cfc 





^k = s. 


'SSk 


- g»a/2 







Sk 




Ck 


- — n 



./2 

ia/2 



= Ckdl + qi 



cos{ipk/2) 

2.3 The integrable motions on S'^{R) 

Since the motion takes place on the sphere, the velocity field is tangent to it: 

Ffc.t = Vfctfe + f/fefifc . 

Futhermore the evolution of the Frenet frame is described by an antisymmetric matrix: 



rfc 
tfc 
life 






XVk 


XUk 




ffc 


XVk 





Ak 




tfe 


XUk 


-Ak 







fife 



or equivalently, by the evolution of the matrix Sk' 



Sk,t - 2 



iXVk Ak + iXUk 

-Ak + iXUk iXVk 



Sk =: TkSk 



(17) 
(18) 
(19) 



(20) 



(21) 



Using Property 2 (i.e. A_t — 0), the compatibility condition 

Edt = dtE , 
where E is the shift operator along the chain: 

Efk — /fe+i , 
applied to the linear problems (|6|) and (pi|), or (|lj) and (|2^), specifies Ak in terms of the velocity: 

A 



Ak 



sin(AA) [l + g2 
A 



^Ffe+i + \-^Uk+i - cos(AA)C/fe 



sin(AA) 



1 + 9^ 
siTupkVk+i + cosipkUk+i - cos{XA)Uk 



(22) 



(23) 



(24) 



(25) 



and yields the kinematics 

2qk^t = -2cos(AA) {Uk+i + qkVk+i) + (1 + ql) [Uk + Uk+2 + qk+l (T4+1 + ^^fe+2)] , (26) 

(1 - ql)Vk+i - (1 + ql)Vk = 2qkUk+i {cosipkVk+i - Vk = sinLpkUk+i) , (27) 

where 

Vk - . ^... Vk , Uk - . ^'..M k . (28) 

sm(AA) sin(AA) 

Using (p7|), equation (pq) can be rewritten in the following convenient form 

2qk^t = -2 COS(AA) {Uk+i + qkVk+i) + n {Uk+i + qkVk+l) (29) 

in terms of the recursion operator 

qk 



7^ := (1 + qt) 



E + E-^ + 2{qk+iE - qk-i){E - ly^ 



^ + <it 



(30) 



which is a real reduction of the recursion operator obtained in |1^, associated with the AL spectral 
problem (p. The operator {E — 1)~^, formal inverse of {E — 1), takes the following explicit form: 

fe-i 
{E - ly'fk = 5] /, + 6 , (31) 

j=ko 

where the choice of the constants fco and b depends on the space of functions fk under consideration. 

Equations ( p^ ( p7| ) show that the kinematics of the curve depends explicitly on the radius of the sphere 
through cos(AA); a drastic way to satisfy Property 5" is to choose the combination {Uk+i + qkVk+i) as 
the real eigenfunction of the recursion operator TZ for the eigenvalue 2cos(AA), i.e.: 

Uk+i+qkVk+i =*fc(A) , (32) 

where 

7^*fe(A) = 2cos(AA)*fc(A) . (33) 

This velocity field, which generates the rigid motion qk,t = of the discrete curve on the sphere of radius 
1/A, can be used to generate more complicated dynamics. For instance, if, on the sphere of radius 1/A, 
one uses the velocity field which would give rise to the rigid motion on the sphere of radius I/Aq ^ 1/A, 
normalized in the following way 

Uk+i + qkVk+i - ,, It^^"^ ,,A, . (34) 

cos(AoA) — cos(AA) 

then one would satisfy Property 3\ obtaining an example of integrable discrete dynamics with sources 
(with a singular dispersion relation): 

gfc,t = *fc(Ao) ^ [7^-2cos(AoA)]gfc,t =0 . (35) 

If, in particular, Aq = 0, then 

A[l — cos(AA)J 



where 



Uk+i-Uk+qk{Vk+i+Vk) = 0, 

Vk+i~Vk-qk{Uk+i+Uk) = 0, (37) 



and 



9fc,t = ^fc(0) = -(l + gD(C^fe+i + C^/c) ^ VKt^Uk+i + Uk 



(38) 



Equations (p^), satisfied by the velocity field which gives rise to the rigid motion on the plane, can be 
integrated to 

Uk = csinilfc , Vfc = -ccosfifc , 



where fifc is the (incomplete) total curvature: 



fc-i 



rji 



y^'-pi 7 (^fc+i 



^t 



"Pk) 



(39) 



(40) 



Consequently equation (|3g) becomes the well-known discrete sine-Gordon (SG) equation [|17[||6|: 

Lpk,t = c(sinrife+i -1- sin fife) ^ (rife+i - ^k),t = c(sinJlfc+i -f sinfJfc) , (41) 

which can also be written in the following form 

(ufc+i - Uk),t = csin(ufe+i -I- Uk) (42) 

through the change of variables 



Ufc-(- 



i-i-ufc = rifc , {uk = yyk-2j) 



(43) 



Using the well-known connection between discrete dynamics and Backlund transformations (BT) |18|, 
we observe that equation (E2) can be viewed as the celebrated BT of the SG equation u^tx — sinu |19|, 
interpreting Uk and Uk+i as two different solutions of it. Therefore equation ( p4| ) has two different geo- 
metric meanings: i) it describes the integrable dynamics of a sequence of points on the sphere, subjected 
to a velocity field which would give rise to a rigid motion on the plane, and ii) it provides an elegant way 
to construct a new (pseudo) spherical surface from a given one. 

The linear dependence on cos(AA) in equation ( p9| ) suggests that Property 5' can be satisfied also 
looking for velocity fields given by suitable Laurent expansions in powers of cos(AA): 



Vk 
Uk 



E (cos(AA))^- 



] = -mi 






, TTll, 7712 > 



(44) 



Substituting equation ( [44| ) into (M) and requiring independence of cos(AA), one obtains a class of inte- 
grable dynamics in the form 

g{'R)qk^t = Hn){l + ql){qk+l " ft-l) , (45) 

where g and h are arbitrary entire functions of their arguments. The simplest examples are the following: 
i) If g = /i = 1 we have 

(46) 



sin(AA) r- „ -, 

Vfc = ^ [tfc - qk-irikl 



A 



qkA = (l + '7fe)(gfe+i -qk-i) 



K 



(1) 



<^ '4^k,t 



, ,fk+l. . fVk-l- 

tan( — - — ) — tan(- 



(47) 



2 ' ' 2 ' 

In this case the velocity field v^ admits the simple geometric construction shown in the figure below, in 
the case of the plane (A — > 0) . 



ii) If g = I and h(x) = x we have 

4sin(AA) r, ,^^, , , - / ,^^, , 1. , 2 ^, ,2 



Vfe = ^ <^ {cos{XA) + qkqk-i)tk+ ( -gfe-iCOs(AA) + ^{l + qk-i){qk - qk-2) - Qk-ilk j fik > , 

(48) 
qk,t = (1 + ql) [(1 + 9^+1)9^+2 - (1 + ql-i)qk-2 + qkiql+i " ^Li)] := i^f • (49) 

As it was observed in llQ], the hnear combination 

qk,t = 4'^ " 2i^f ^ (50) 

of the first two flows reduces, in the continuous hmit, to the mKdV equation 



iii) If g{x) — X and h — 0, then 



tan(AA) 
Vfc = ^ 



^k ^k + Uf. rifc 



and the tangent velocity 



(52) 



qKt^-\{i + qi){uU-\>~uir\>) , (53) 

where the normal velocity field satisfies the difference equation 

qk.. + -^) ui^ + (,. - 1) ui^ - L,, - -^) ut'^ - (,. + 1) f/<:;) = , (54) 

qk+ij \ qkj V 9fc+i/ V qkj 



2^r^--(.. + ^)c/i;P + ('z.-^)t/l.:l^ (55) 



We remark that this integrable dynamics corresponds to the velocity field 

C/fc+i + qkVk+i = jT-TT : (56) 

cos(AA) 

which would give rise to a rigid motion on the sphere of radius I/Aq — 2A/7r. 

3 The motion of a discrete curve in S^{R) 

3.1 The discrete curve in S^{R) 

The plane tt^. is, as in the previous Section, the oriented plane given by the ordered pair of vectors 
(Ffc , ffc+i). Similarly, the unit vectors tk and tj, are tangent to the big circle passing through the points 
ffc and r^+i. 

The tangent plane Tk in the point fk to the discrete curve is defined by the two vectors i'k-i and 
tfc. The angle of geodesic curvature ipk is the angle between the vectors t'^, and t^+i. The orientation of 
the plane Tk+i is choosen by the condition < (pk < n. The normal vector n^ is the unit vector of Tk 
orthogonal to t^ and such that the pair (t^, fifc) gives correct orientation of the plane t^. 

The binormal vector b^ = fj, x tj, x nj,, together with the normal vector fifc, define the normal plane 
I'k in point r^. The orientation of the plane I'k is given by the ordering (nfc,bfc). We first remark that 



the vectors f^+i and t'^ span the plane tt^; moreover, since tj, e Tfe+i, it follows that bfc+i±t'^ which, 
together with the obvious condition bfc+i_Lffe+i, imply that b^j+i E i^k = ^k ■ The angle 9k between the 
vectors b^ and bfe-|_i (which belong to the same plane i^k), measured according to the orientation of the 
plane i^k, is called the angle of (the geodesic) torsion of the curve in point ffe. In the case of R"^ it is 
ilustrated on the figure below 



In the continuous limit, the angle of torsion 9k reduces to the familiar torsion t{s) of the corresponding 
continuous curve in the following way: 

— — > t{s) . (57) 

3.2 The Hasimoto transformation and the spectral problem 

As in the case of the continuous curve, it is convenient to change the basis of the normal plane Vk through 
a rotation of the angle tJfc-i = {E — l)^^9k ~ J2i ^ ^i- 

fifc = coscrfe_i fifc - sino-fe_i bfe , (58) 

h\ = sinak-i n^ + coscrfe_i bfe . 

This change of basis corresponds to a partial "integration" of the Frenet equations in the normal plane, 
since the vectors fi^ , n| do not vary from the point of view of the normal plane (this is the discrete analog 
of the parallel transport in the normal bundle) . 

As in the case of a continuous curve, it is convenient to interprete any vector of the normal plane Vk 
as the complex number 

$k = Re0fe n^ + Im^fe nl ^ (jy^ ^ Rc(/>fe + i Im^^ • (59) 

The transition from the basis (f^, t^, n^, n\) to the basis (ffe+i, t^+i, n^^-^, ii^k+i) '-^^ ^^ obtained by the 
superposition of the following five rotations: 

i) rotation in the plane Vk of the angle dk-i 

(ffc,tfc,n^,n'fe) — > (ffc,tfc, fife, bfe) 
ii) rotation in the plane TTfe of the angle 2a = AA 

(ffejtfe, fife, bfe) — > (ffe_|_i,tfe, fife, bfe) 
iii) rotation in the plane Vk of the angle 9k 

(ffe+i,tfe, fife, bfe) — > (ffe+i,tfe,fife,bfe+i) 
iv) rotation in the plane Tk+i of the angle (pk 

(ffe+ijtfe,fife,bfe_|_i) — > (ffe_|_i,tfe+i,fife_|_i,bfe+i) 
v) rotation in the plane Vk+i of the angle — crfe 

(rfc+l,tfe+l,fife+l,bfe+l) > (ffe+l,tfe_|_l,fife^;^,fife_|_;^) . 



Rotations in R^ can be represented in a similar form as in the previous Section. A fixed Cartesian basis 
in R^ can be identified with the following 2x2 matrices £ = {eo = Id, ei = ias, £2 = ici, £3 = —102)- 



vGR^ ^ (t(v)=^v'£, ; ||v||2=dcta(v) 



(60) 



1=0 



Any rotation O can be represented in terms of two matrices P and S of the SU{2) group in the following 
way 

w == Ov 4=> cr(w) = P" V(v)S' . (61) 

When we are given a sequence of orthonornial frames J-'i^ = (e^, fj., gj., hj.) and the corresponding matrices 

Pk, Sk 

Tk = Pk^SSk (62) 

then we define matrices Ak and Bk by the relations Sk+i = AkSk and Pk+i — BkPk] moreover 

PkTk+iS^' ^ B^^SAk . (63) 

If the rotation takes place in the space perpendicular to the first vector of the basis J-'k (i.e. e^+i = e^) 
then Ak — Bk- For example the matrices 



^-^e/2 

C*«/2 



Ak — Bk — 

represent the rotation in the plane (gfc, h^) of the angle 9, and the matrices 

Ak = Bk = 



cos{ip/2) sm{ip/2) 
— sin((p/2) coa{(p/2) 



(64) 



(65) 



represent the rotation in the plane {ik,gk) of the angle (p. 

We will also need the representation of the rotation in the plane (e^jf^) of the angle 2a, given by 

matrices 

" e"*" 



Ak 



e*" 
e'" 



, Bk 







(66) 



We are now in the position to calculate the transition matrices when the sequence of the orthonornial 
frames is the sequence of the frames along the discrete curve Tk = (ffc, t^, n^, n^. The corresponding 
spin matrices Pk, Sk are subjected to the relations 



Pk+i = 



C°'cos{ipk/2) e"^sin{ipk/2y'"' 

-e-'" sin(^fc/2)e-^'"'» e*" cos((/?fc/2) 



Pk 



1 



v/T+kP 



(-' (qk 
-C'qk C 



Pk =■■ BkPk 



(67) 



where 



Sk+i 



e'" cos((pfc/2) 
-e*" sm{ipk/2)e~ 



1 



VT+W' 



C C'qk 
-Qqk C 



^■sin{ipk/2)e""' ' 
*"cos((^fe/2) 

^k ='■ AkOk , 



qk=taniipk/2)e"'^ , C = e'' 



Sk 



(68) 
(69) 



10 



We remark that the transformation (p9) between the curvature and torsion angles and the complex field 
Qk is the discrete analogue of the Hasimoto transformation, introduced for continuous curves in [ pO[ . We 
also remark that Bk = ^fcdiag[C~^, C^]- 

As before the linear problem can be transformed into the AL spectral problem: 



Sk+i 



C qk 

-qk C 



Sk , Pk+1 



C qk 

-qk C 



Pk 



(70) 



(but now in the reduction r^ — —qk) relaxing the conditions detSfc = detPfc = 1 (but keeping still 
detSfe = deti-fe) and changing the basis (ffe,tfe) of the plane TTfc to the new basis (f'^,afe). 

3.3 The integrable motions in S^{R) 
Proceeding as in Section |2.3| we have that 

sin(AA) 



rfe.t 



( Vfetfc + (f>k] , 



(71) 



= Re^k nl + Im^k n^ = Rc{<j>ke-"''-') hk + lm{<j>ke-"''-')hk 



and, consequently. 



Sk.t — TkSk 



(72) 



Tfc = - 



e su{2) 



Ik Sk 
4 -7fe 

The compatibility condition between equations (|68|)([7l|) and ( |72| ) specifies the entries 7^ and Sk in terms 
of the velocity field: 



Sk = < ^'/'fe - «(0fc+i + qk{Vk + Vfe+i)) 
Ik - Wk+sm{XA)Vk 
Wk+i~Wk = ~Re{iqk[(l)k+2-(t>k+qk+i{Vk+i+Vk+2)]) 

and yields the kinematics 

2qk,t = -2cos(AA) {4>k+i + qkVk+i) + (1 + l^fel^) [0fc + (f>k+2 + qk+i{Vk+i + Vk+2)] 
+ 2iqkWk + qk{qk4>k - qkipk) 
= -2cos(AA)(0fe+i +(7fcVfc+i) +^((/'fc+i +g/cVfe+i) , 

(1 - \ql\Wk+i - (1 + kfe^Vfe = qk<i>k+i + qk<j)k+i , 

where 

njk := (1 + IgfcH f/fc+i + fk^i + 2{qk+iE - qk-i){E - l)-^Re '^f" ) 
- 2iqk{E -iyhm{qk+ifk -qkfk+i) ■ 
Substituting the ansatz 



Vk 

4>k 



^(cos(AA))^ 
3=0 



V, 



(m~j) 



li^m-j) 



(73) 



(74) 
(75) 



(76) 
(77) 
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into equation ( [74[ ) and requiring independence of cos(AA), we finally obtain the following class of inte- 
grable dynamics: 

qk,t ^ ho{TZ){l + \qkf){qk+i - Qk-l) + hi{n){iqk) , (78) 

where ho and hi are arbitrary entire functions with real coefficients. 

We remark that, in the degenerate case of the curve on S'^ {qk G R) we are forced to choose hi — 
and only the first hierarchy survives. 

We also remark that equation d7J) implies the following interesting connection: 

Kt'-^tl+l'^y^Tl (79) 

between the integrable commuting flows 

^(m) ^ ^,n-i(^ ^ \qkf)iqk+i - qk-i) and/or K^^ = 7^™(igfe) , m > (80) 

and the velocity fields. In the continuous limit this reduces to the result of | |2l[ |. 

The simplest examples are the following: 
i) liho = 1, hi = 0, then 

sin(AA) ,- ^ , sin(AA) ,- , ,^ ^ ,„,n 

Vfc = — ^(*fe - ^k-i) = — ^(tfc - \<lk-i\&k) , (81) 

qk,t ^ (l + \qk\^)iqk+l - qk-l) =■■ Ki'^ ■ (82) 

ii) If ho{x) — 1', hi = 0, then 

4sin(AA) ( f , ,, I, _ _ xA - - 1 

Vfc = <^ I cos(AA) + -(qkqk-i + qkqk-i) \tk + q>k} , (83) 

(j^k = ( -9fc-i cos(AA) + 2 [(1 + \qk-i?){qk -qk-2) - qk-i{qkqk-i +qkqk-i)\ 

qk,t = 4'^ := (84) 

(1 + \qk?) [(1 + \qk+i?)qk+2 - (i + \qk-i?)qk-2 + qk{ql+i -ql-i) + qk{qk+iqk-i -qk-iqk+i)] ■ 

The linear combination 

qk,t = 4'^ - 2^'^ (85) 

reduces, in the continuous limit, to the complex mKdV equation 

V-,*' = i^,sss + lm^i^,s , i^is) :- hm ^ - K{s)e'^^'^ , t' = 2A3t . (86) 

iii) If /iQ = and hi(x) = x, wc obtain 

^ 2sin(AA) ^.^ 2shi(AA) .(^fc-i.r .„„. 

Vfc = ^^ (iQfe-i) = 3;^ tan(^— )bfc , (87) 



<7fc,t = «(1 + |gfcn(9fe+i + 9fc_i) . 
The following combination of equation ( |88| ) with the "zero order flow" gfc j = iq^ ■ 

qk,t = i [qk+i - "^qk + qk-i + \qk\'^{qk+i + qk-i)] (89) 
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reduces llG], in the continuous limit, to the NLS equation 

i^.t = ^.ss + ^\^\^tl^ , t' = -Ah . (90) 

which describes the motion of a vortex filament in the localized induction approximation ||22l[Q. We 
remark that, in this approximation, the velocity field which governs the motion of the vortex depends on 
its curvature k through the relation 

r t = Kh ; (91) 

therefore since equation ( p9| ) has in R^ a velocity field of the same type: 

rfc,t = 2Atan(^)bfc , (92) 

we expect it to be a good candidate for describing the motion of a discrete vortex in the same approxi- 
mation. 
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